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ABSTRACT 


Simultaneous measurements of wave height and two orthog- 
onal water particle velocities were made at six elevations 
in nineteen meters of water using a penetrating wave staff 
and an electromagnetic flowmeter. Moderate swell and low 
meee CONGItTIOons prevailed during the Sen iriene The mea- 
sured wave-induced velocities were two to four percent 
Breater than those predicted using linear wave theory. 
Coherence of the wave height and wave-induced velocities in 
ie ssteniticany energy-density range was computed to be over 
0.85, indicating that the motion was almost totally wave 
meemced. Av higher frequencies iC was apparent that the 
motion was primarily turbulence. Phase spectra computed for 
the measured wave heights and orbital velocities compared 
very well with linear theory. Measured frequency distribu- 
tions were compared to both Gaussian and Gram-Charlier 
distributions using the chi-squared goodness-of-fit test. 
Qualitatively, the Gram-Charlier distribution gave the 


better fit to the data. 
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i Se RODUCT ION 


A. BACKGROUND 

Measuring instantaneous water particle velocities in 
Meemoresence of wave-induced motion has been difficult due, 
permearily, to a lack of instrumentation. This is particu- 
larly true for measurements made in the field. Until re- 
Centvily, there has been no adequate instrument available to 
measure rapidly fluctuating flows found beneath wind waves 
and swell. As well as the usual requirements of reliability 
ana accuracy, such an instrument must also be rugged, have a 
low response time, and have good directional measuring 
emaractleristics. 

Earlier attempts at correlating water height and particle 
motion have generally been accomplished in shallow water next 
Memmne  GMOre, OLl pliers, or in estuaries. Historically, 
Inman and Nasu (1956) made one of the earliest attempts at 
field measurements of wave-induced water particle motion, 
using a force meter in shallow water and compared their 
results to solitary wave theory. Miller ane Zeigler (1964) 
used both an acoustic Rete (aaa ONneune GCoOppler shane 
principle) and an electromagnetic flowmeter to measure water 
weet icle movion in the surf zone. They compared their 
measurements to higher order wave theory and obtained 
qualitative agreement. 

Nagata (1963, 1964a, 1964b) used an electromagnetic 


flowmeter ina series of experiments to measure water 


EG 





@emunmecle velocities in the surf zone and to determine wave 
directional spectra. Bowden and White (1966) and Simpson 
(1968) performed field measurements using an electromagnetic 
Plowmeter developed at the National Institute of Oceanography. 
Like Nagata, Bowden and White (1966) used water particle 
velocities to determine the directional spectrum of waves. 

Simpson's experiment was performed at an open coastal 
site in water approximately 6 meters deep. The waves, 
w@mmmoar. ly wind penerated sea waves, were developed over a 
mepeine fetch of cl0 km. Simpson employed spectral methods 
in the analysis of the waves and of the wave-induced water 
particle velocities. His results were consistent with 
Hinear wave theory. 

Seitz (1971) carried out a similar experiment in that he 
measured sea surface displacement and water velocity compo- 
Me@es an an eStuary at a depth of one meter. He compared 
the measured velocity spectra to a derived velocity spectra 
Gevained from direct measurement of the surface wave field 
Using linear theory. Seitz had to contend with a shift in 
moe observed frequencies, the napaedas of which was a 
function of the tide induced velocity. Since at the fre- 
quency of peak amplitude the error introduced was approxi- 
mately 6.7%, he did not correct the data. The total water 
particle motion at depth can be described as the sum of The 
turbulent velocities and the surface wave induced motion. 
Seitz attempted to separate the two based on the slope of 


the spectrum in the inertial sub-range. 


dha 








Mhe lack of instrumentation has not been the only 
hindrance to wave height and velocity studies. Most research 
has been conducted at shallow water inshore sites because of 
thie general unavailability of a deep water, exposed instru- 
ment platform. Shonting (1966) used the Buzzard's Bay 
Peerence Light Station on the South Coast of Massachusetts 
Memoovain a Cwenty meter Water depth at an open ocean site. 
Pewused ducted current meters to measure the orthogonal 
Semeonents of the water velocity under waves. His experiment 
was designed primarily to evaluate these meters. An inherent 
Seer acuity with any type of ducted meter is their poor 
directional measuring characteristics when the angle between 
Meecr axis and the flow is large. 

Bordy (1971) made measurements in 20 meters of water 
from the Naval Undersea Research and Development Center (NUC) 
Oceanographic Research Tower located approximately one mile 
off Mission Beach, San Diego. He compared wave height and 
water particle velocity using spectral analysis, and found 
a general agreement with linear wave theory except for 
measurements of phase spectra at ane deepest measurements of 
Gen meters. The phase shift encountered was ule =p faa san 
the phase angle error analysis. During the experiment, the 
NUC tower's permanently installed wave staff system was used 
for water level sensing. The flowmeter was removed horizon- 
tally from the wave staff by approximately five feet. It is 


felt that due to this displacement, an error was introduced. 


re 








During 8 June 1972 records of orthogonal water particle 
velocity and instantaneous water elevation were again 
measured at the NUC tower. The data were taken in conjunc- 
tion with an experiment involving the investigation of wave 
meemeed fluctuations in acoustic phase shifts, temperature, 
emoesalinity. 

The results of the wave and water particle motion 
Weeourements are described here. An electromagnetic flow- 
meter was used to gather the wave particle velocity data 
at several depths. These measured values of wave-induced 
motion were compared with values predicted using linear 


weer y and the wave height record. 
Bees obARCH OBJECTIVE 


mae purpose of this investigation was four fold: 

1. To measure Simultaneously wave height and two 
Gmmmoronal wave-induced water particle velocities. 

flO Pperiorm a spectral analysis on the data to obtain 
imeem Svavistical properties. 

3. To spectrally compare the data to first order wave 
Theory. | 


4, To determine the statistical geometry of both the 


sea surface and the measured velocities. 
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TT. EXPERIMENTAL SITE AND INSTRUMENTATION 


Peet XPERIMENTAL SITE 

The experiment was conducted at the Naval Undersea 
Research and Development Center Oceanographic Research Tower 
at San Diego, California, 6 and 9 June 1972. The NUC tower 
meme ocalred approximately one mile off Mission Beach, San 
Diego in a mean water depth of 19 meters. It is constructed 
of steel and concrete and is permanently imbedded in the sea 
floor, Figure (1). The forty foot square tower has three 
i7ers, The first is used primarily for gaining access, 
while the second houses machinery and provides work space 
ae lacing instrumenteewim thie water. The third, which is 
eeemesed, 1S tised for tUhewmeuciamea ve: electrical equipment 
meemotso provides sleeping space tor six men. The legs of 
meemuUOwer Sitepe Tive decrees wa vair of tracks is mounted 
Oemeie North, Seu, and westssides fer lowering scientific 
instrument packages into the sea to any depth Aneel cine 
Secametloor. The western side of thewtower was used for the 


experiment because of its exposure to the predominant swell. 


B. INSTRUMENTATION 

Water particle velocities were measured with an Engineer- 
ing Physics Company water current meter Model EMCM~3B. The 
sea surface elevation was measured with a Baylor Company 
wave staff system Model 13528R. The penetrating wave mea- 
suring system gives a direct measure of the sea surface 


elevation. 
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Figure 1. Naval Undersea Research and Development Center 
Oceanographic Research Tower, San Diego, California 
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The model 13528R Baylor Company Wave Staff System is an 
instantaneous water level measuring device consisting of 
two tensioned 1/2", 6x19 IWRC 304 SS wire ropes. Together 
with these lengths of stainless steel wire rope, the trans- 
mececr produces an electrically linear direct current output 
that is proportional to the amount of wave staff above a 
seer, Circuit produced by the water surface. System aecuracy 
is stated as within one percent of actual height and the time 
constant is less than 0.06 seconds. 

mo .srecei ved Trom the mManimacuurers., The wor. ine = bomen 
@, tne staff was 50 feet; however, in order to take water 
barcicile velocity measurements near the surface direcvily 
macder the staff, 10 was shortened to 20 feet. After 
emeorcvening, the staff was recalibrated. The resus 
calibration curve is shown in Figure (2). The rails on 
which the instrument cart rides slopes with the Tegswei tie 
tower. This necessitated the rigging of an aluminum "I" 
beam with a trolley to support the Baylor wave staff system 
directly over the electromagnetic flowmeter. The wave staff 
was tensioned eo sadionaie the weight of the INStPUMeENL "Care 
and frame to be exerted Bi Ene lower end of the wave staff 
cables. The estimated tension was approximately 800 pounds. 
The instrument and installation are shown in Figure (3). | 

The EMCM-3B electromagnetic current measuring device 
measures two orthogonal components of water particle velocity 
through a range of zero to five M/SEC. Other characteristics 


of the flow meter are: 
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Figure 3. Model 13528R Baylor Company Wave Staff system 
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Output Voltage: +5 volts full scale 


Electrical Time Constant: 0.2 seconds or longer, 
adjustable by user 


Breror Band: the largest of 
1) +1 percent of full scale 
2) +5 MM/SEC 


L 
3) 5/(T*) MM/SEC 
where T is the time constant 


Transducer Size: 283740. = 26 as 


mide Operation of the unit is based on Faraday's principle of 
magnetic induction. A current is induced when a conducting 
fluid flows around the probe through an alternating magnetic 
field generated by a toroidal magnet inside the instrument. 


The equation describing the production of current is: 
B= u XB, 


where EF is the induced voltage, u the Vel OC iy 7 -Ol une 
conducting medium, and B the magneric field.) The snduced 
ereectromovive force is linearly proportional to the amount 

of fluid flowing through the reene ce ttenee The magnetic 
fone ld, B, is most intense in Oe aac vicinity OL the probe 
and decreases according to the inverse square of the distance 
meom the instrument. Accordingly, the flow velocity near 

the probe is weighted most heavily in its contribution to the 
measured velocity. Qualitatively, only that water which is 
Within a distance of two to three probe radii contributes 


SEchaticanviy to the flow induced volvage. 


iPS, 








mie tlow meter was originally calibrated under steady 
flow conditions by the manufacturer. Different response 
Characteristics are to be expected under unsteady flow 
Conditions; consequently, the probe was recalibrated by 
oscillating it in a water tank. The equipment schematic is 
shown in Figure (4). The carriage on which the flow meter 
probe was mounted travels back and forth on rails. It is 
Griven by a variable speed electrical motor. The peak 
@ereriase velocity was calculated from the tangential 
ipeocicvy Of the motor arm. The ratio of carriage velocity, 


V meme Velocity measured bY The insprumens., tas was 


ee 
meuna ior different angular velocities. This was done for 
each pair of electrodes by orienting them parallel to the 
mwow 17) turn. 

The instrument's response was very linear with respect 
memirequency up to at least 0.5 Hz. Beyond this frequency 
The cart upon which the probe was mounted became unstable. 
Wisc. 5 Hz was the maximum frequency investigated during 
Calibration. This limitation was not deemed a debilitating 
mrereor Since 0.5 Hz was above the upper frequency limit of 
Significant wave energy. The results of the calibration are 


shown in Figure (5) and the instrument and installation in 


Figure (6). 
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Figure 6. Engineering Physics Company Model EMCM-3B 
Electromagnetic Flowmeter Installed in Position (arrows) 
Gr the InsSerument cart 
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sea surface elevation and two orthogonal components of 
Wower particle velocity were measured simultaneously. The 
Cart on which the flowmeter was attached was raised and 
lowered so that measurements could be made at various depths. 

Data irom the flowmever and wave gauge were recorded on 
a Sangamo Model 3500 fourteen channel FM tape recorder 
GObring the ten data collection periods. The data chosen for 
Peelysis was Obtained during the first six recording periods 
Bemeo June 1972. Each processed run was about 20 minvtes 
m@eeacion. After initial processing, the other four dara 
@ollection periods were deemed too Short to be statistically 
representative. 

Summarized in Table I below are the velocity components, 
Water depth, depth of the instrument, and record Lenieune et ore 
eee Of the processed runs. Instrument depth refers to the 
Gepth of the flowmeter below mean sea level. The u-component 
oF velocity was oriented north-south approximately parallel 


moerme shore, while the v-component was measured east-west. 
Cerrect flowmeter srLonenet on was obtained using a plumb-bob 
and level. The wave gauge was always maintained directly 
over the flowmeter. As the cart was lowered the top of the 
wave staff was pushed further away from the tower to compen- 


sate for the five degree tilt of the tracks to which the 


cart was fastened. 


a4 








mime Wave Stari's perpendicularity to the horizontal was also 
ascertained by the use of a level. The wave staff installa- 


Mion is pictured in Figures (7) and (8). 


Ste jaa b 
DEP DEPTH OF COMPONENT LENGTH OF 
RUN (meters) INSTRUMENT OF RECORDS 
(meters) VELOCITY (minutes) 
il io eal 7.04 Dacha, 24.96 
2 ores 2 T. 62 u & Ww Ceol 
S HO5.25 eon u& Ww el. 35 
ones 10.69 u & wW 24.85 
5 19.32 8.66 u & w 18.88 
6 19. 36 | Sw u & w 19.52 


HWhe sea surrace was in general disturbed by swell coming 
from the west-southwest. The sky was overcast with a low 
Gemauvus layer and little wind, less than 10 knots generally. 
wea surface temperature was obtained using a bucket thermom- 


eter. Table II summarizes the environmental data. 


TABLE II 
TIME OF AIR SEA WIND CLOUD SWELL 

RUN START OF TEMP SURF DIR Vin CO VERO Db 
RUN (PDT) (°C) TEMP (°C) (°T) (knots) (%) 

iL 1405 IR A alee axe 20 8 1¢c0 36s WS 

2 1600 19.6 18.95 250 10 100 WSW 

3 1654 epee Ts 6 20 9 100 wWSswW 
1800 Te. 7 a bono 225 6 100 WSW 

5 1830 Tee) Shesos 225 7 100 WSW 

6 aoweye, 7 Oar 7 235 \ 100 WSW 


25 











Figure 7. Model 13528R Baylor Wave Staff Installation 
Showing "I" Beam and Trolley for the Top of the System. 





Figure 8. Pulley Attachment Point for the Bottom of the 
Baylor Wave Gauge (arrow). 
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Table III summarizes the temperature profile data taken 
eueeme beginning of each data collection period. This data, 
which shows a remarkable variation in layer depth and 
intensity of temperature gradient, was obtained using a 
Meenmanical BT. The NUC tower was originally constructed to 
study internal waves which are prevalent in the area. The 
variability in temperature profiles is attributed to internal 


waves. A typical temperature profile is shown in Figure (9). 


a 


TABLE EEE 

Br TOP SOr BO Onor Or OF BOTTOM 
men 6SURE 18t LAYER 18 LAYER ena TAYER TEMP 

EMP (oC) Ce7m) ey) (267) (Se) 
a 2280 Pees val Ore ele) ex x 15255 
2 22.5 222/30. Wong ee CO 27 15.6 Sem lomo 
3 E53 23.3/6.00 “2G 7s. POURS 16.46 ile 8) 
y ee. OTe 1/4. Oy Ome Gr aa: none ley] 66 16.66 
5 Bec Pl.6/ 5.409507 FORGO, mens eee Lon 33 
6 ag, 9 19:,9/5.79 16217 15.24 Sex Gi, ae 


jie analog FM recordings were digitized using a hybrid 
system consisting of a Scientific Data Systems (SDS) Model 
Se00 analog computer and a Xerox Data Systems Model 9300 
(uescal computer. For comparison purposes the data were 
transcribed onto a rectilinear eight channel Clevite strip 
chart recorder after a tenfold amplification through the 
(SDS) 5000 analog computer. An example of the transcribed 


Geta cerberedio1gizareom, aS Saown in “igure (10).. Yageer 
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Figure 9. Typical Temperature-Depth Profile. 
Run Number 4. 8 June 1972. 
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res rizing onto the seven track tape used by the 9300 system, 
the data were transcribed onto a nine track tape using the 
Naval Postgraduate School's IBM Model 360 digital computer. 


ims ecOMPDUutTer was also used for the statistical analysis. 
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A. LINEAR WAVE THEORY 
Mae, clevation Of tune Surtace Nit) can be described as 


the superposition of an infinite number of sinusoids of the 


ierein ¢ 
co re . CO 
n(t) = ) a, cos Ogre -c. item) = ) ni, Ciesae) 
n=1 n=1 
where x is the horizontal Cartesian coordinate, t is the 
iALme i HS) a NOrizontal vector wave number, 0. the trequene, 


Y) 


mebeaced in linear theory to k Diy: 


Cun : 
oes 30 ie volg ic (yn) 


where g is the acceleration of gravity, a, is the amplitude 
and E, the phase angle. n(t) represents the sum total of 
ali component wavelets. Summing in the manner of Equation 
@r)) implies a linear system and restricts the analysis to 
the use of linear or Airy wave theory to describe the wave 
Haowced mOouron . 

Linear wave theory can be used to Deed Cl. Tae Wave 
induced water particle velocities. The equations for the 
horizontal and vertical velocities respectively are: 


aek, cosh k , (atz) 


> > 
we cos(k,+x-o,, tte) 


T 
i138 


1 6. Sinn eon 
n n 
gk cosh k,, (htz) (4.3) 


or = y ee. Seeeeenes 2 eines nN 
=F oF Sinh Kf n 


eel 


and 


ach alen Saliglen ) cml lonrs | 
ad Yaa n Py aaa 
cu) = ) aaa aa oc ee sin(k, HOG) oe) 
n=l n n 
(4.4) 
_ ; gk sinh k,, (ntz) ; 
= Oo sinh k_h n 
n=l n 1g eE=90° 
meeweo"h™“is“the totelewaver depth and zis the depth of 
gaverest measured positively upward from the still water 
tevel. 
The terms 
gk cosh k_(htz) 
iso Sola ae. (4.5) 
oj Salata eee’ 
n n 
and 
gk sinh k_ (htz) 
is = (4.6) 


ee ee ne Ee nee 
Cian lel aia 3 
n n 


jMemeescent the transfer functions relating the horizontal and 


vertical velocities to n(t). 


B. LIMITATIONS OF FIRST ORDER WAVE THEORY 

In the formulation of linear wave theory, the boundary 
conditions are linearized in order to obtain an analytical 
souucion. in the linearization, it is assumed that the 
emplitude, a, is small compared to the wave length, L, that 
mona /l << 1, 

Higher order solutions to the boundary value problem, 


generally obtained by perturbation analysis, give a better 


oe 





representation of a constant profile wave. However, the 
nonlinearities Haire dsad ie como OU rhOnmoree Ludes tile LTr 
Maeswnere the principle of superposition is invoked. This 
im turn eliminates their generalization to a stochastic 
process. 

Dean (1968) analyzed the various wave theories with 
mespect to how well they fitted both the dynamic and kine- 
matic free surface boundary conditions. He showed that 
iment Order theory provided a good solution to the wave 
equation for waves of permanent form having small steepness 
in intermediate water depths (the assumption of linear 
theory stated above). 

mime condivions prevailing during this experiment fella 
within these limitations. The rms wave height averaged 0.27 
meters with a mean period of 16 seconds corresponding to a 
deepwater wavelength of 209 meters in 19 meters of water. 
Hence, if lineer theory provides a reasonable approximation 
for the average wave conditions, the superposition of a 
large number of components, each described by J3near Taecry. 


mieuLd provide an even better description. 


! 
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V. THEORETICAL CONSIDERATIONS OF DATA PROCESSING 


BY. SrATISTICAL ANALYSIS 


The digitized records were spectrally analyzed using 
the methods of Blackman and Tukey (1958) and utilizing the 


IBM Model 360/67 computer at the Naval Postgraduate School. 


fee COMPUTATION OF THE ENERGY-DENSITY SPECTRUM 

The covariance function ets) of time series randozn 
data describes the general dependence of the values of the 
Gaea av one time on the values at another time for the same 


eee Set. in integral form 
lim 1 
me (1). = Se SG eae) or Sil) 
: 0 


which is the average of the product of the value of a data 

mount o(t) with that of an adjacent point at timest later. 

The resulting average product will approach the exact auto- 
@errelation function as T approaches infinity. 

When using a covariance function generated with a finite 
@aue set, it is necessary to apply a lagw@window To oT). 
Peecarzen window, as used “in this analysis, Nas the desirable 
preperty of no negative side lobes thus eliminating any 
moscoitble numerical instability which mijene arise=due To 
resonance effects of the side lobes. This is particularly 
critical in computing cross-spectra. The Parzen window is 


defined as: 
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P(t) = 1 - 6(1/T_)¢ + oye eee rae eee 
= 3 
ce C/A) eee ee 
= 0 7 eee 
m 
Coe) 


where Dn eo he eval dec mame . 
ine energy-density.specurmum as cCempuLred from Eneerours er 
imeamstorm of the covariance function with the Parzen window 


meee d to correct for finite record length is 


ee = {| P(t) (7) ae oe : (S25) 


CO 


were f denotes a particular frequency. 


ie energy-density spectrum is the statistical estimate of 
Miemwenergy-density distribution in the Various frequency 
bands comprising the spectrum and the area under an energy 
Belisaty Spectrum is equal to twice the variance of the random 
variable. With respect to the wave height spectrum,the area 
tomproportional to the total potential energy-—density of the 
sea surface elevation, while the area under the velocity 


spectrum is proportional to the total kinetic energy-density. 


Bree COMPUTATION OF LAE CROSS=SERCTRAG DEVvSuly 
The cross-covariance function is similar to the covariance 
fucntion with the exception that it compares two time records 


laeeed Dy a tCime interval 1. -Thus, an 2ntvecgrad- form 


do  (t) = es J Pat) acc et) der Cio) 


ji x and y are statistically independent and the mean of 


Maem x and y is zero, ee, =O atOreuc le! aes ener. Peet) 
memequal to the product of the means. The lag window 
wemoLragerations for a finite record length are identical to 
wese for the covariance function. The cross—spectral 
density is therefore the Fourier transform of the Parzen 
Window lagged cross-covariance function 
a i. : -ieTmft 

ett) = Je ay ee) e dt Co) 

ine Cross=—spectrum can be defined in terms of its real 
weeerimaginary parts Since, unlike the energy—density spectrum, 


ieees Not an even function. Thus 
o(f) =O, (£) - 1. (2) (5.6) 
where the real part of 2 y(t) is the co-spectrum 
© 
Cy (f) = 2 J Co,y(t) + o,.(-1)] cos(2nft) dt , (5.7) 
and the imaginary part the quadrature spectrum 


0 


OA ea 2 J Co, ,(1) - o,,(-t)] sin(2nft) dr. (5.8) 
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fee COMPUTATION OF COHERENCY 
A useful quantity arising out of the computation of the 
@ness— and energy-density spectra is the coherency. 


Coherency is defined as 


[o. (f) 1° 
Vy (f) ~ a 7 (Dy 
Oe (f) O(L) 
where 
Bey (f) <2. (5.10) 


When AR = 1, for ali-f£;, thestwo records x< end. 7 farce 
fully coherent and statistical dependence can be assumed. 

On the other hand, if Vy (Lf) = 0 for all f, x and y are 
statistically independent. Bendat and Piersol (1966) 
mmeeibputve coherency valves Less than unity tc three possipee 
causes: the system is not linear, Sere eos noise is 
present in the measurements, or the output response is due 


wommore than one input function. 


D. COMPUTATION OF PHASE ANGLE 

The cross-spectral phase angle represents the average 
angular difference by which the cross-correlated components 
of y(t) lead those of x(t) in each spectral frequency band. 
The phase angle can be computed using the co- and quadrature- 


Becerura 


Exy \f) tarot LQ, ff) vi ee : (5.11) 
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fee ROBABILITY DENSITY FUNCTION OF THE SEA SURFACE 
AND INDUCED WATER PARTICLE FLUCTUATIONS 


iiee probability density function, pix). of randomecdaccs 
describes the probability that the data will assume a value 
within some defined range at any instant of time. 


Mathematically 


BY cry ee ae) «eo 
Where T, is the amount of time that ¥(t) fale ine vie 
range (x, xtdx) during an observation time T. 

for Chis experiment, the» sea surtace isp lacememyaiae 
assumed to be the sum of an infinite number of independent 
memeonents of random phase. Under these conditions the 
Measured frequency distribution of the sea surface displace— 
meme CLeENGS to be Gaussians lHoeme zero mean the Gaussian 


meeeapl. lity density function 18S gBiven by 


ea/e 


po(—é) = = (5.13) 


Cu 


where € is equal to the value of the measured quantity 
divided by the standard deviation. However, if the non- 
linearities of the boundary condition become important, 
departures from Gaussian are expected. Longuet-Higgins 
(1963) showed for "moderately" non-linear swell in deep 
feacer thas the provebllity WistribuciCm Of Weis eurer 
described by the Gram-Charlier probability density function. 
Pec (é) is given by: 
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3 My, ma” 
Peo (6) = pal) [1.0 4+ sy Hy + Gop Hy + ae He) + ...J. (5.14) 


Here the m.'S are cumulants with ue being the skewness of 
the data record, mM), the kurtosis minus three, and oe the 
Hermite polynomials of degree n. The Hermite polynomials 


ere generated from 


yy eh = Bape) EE nna yfncadineg) SPA a5) 
Teeforé 

Hp = 6° -3 , 

i = e" ~ 6E° eons 
a ba 6 ~ 15e! : Naee - 15, (5a) 


iMeeche distribution of € is Gaussian, the skewness is zere 
and the kurtosis is equal to three which results in these 


and all higher order cumulants being Zero making 


Pag (6) = pglé) 


Kinsman (1960) made measurements of wind generated waves 
and showed their distribution was only approximately Gaussian 
and better approximated by the Gram-Charlier distribution. 

He found the distribution exhibited a positive skewness and 
attributed this skewness to the physical fact that waves have 
@ Gisvince. tendency ve, torm short. Sharp crests andegicore 


Shallow troughs. 
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FF, GOODNESS-OF-FIT TEST 

The chi-square goodness-of-fit test is used to compare 
a measured probability density function (PDF) to a theoret- 
meal one. A statistic with an approximate chi-squared 
distribution is calculated as a measure of the discrepancy 
between an observed PDF and the theoretical PDF. The 


chi-squared parameter is: 
i 
Sa Soe 


where r. 4S the observed fregucncy ol voccurence . BR, ete 
Maeerelical frequency of occurence, ana J is the number om 
intervals into which the frequency distribution is divided. 
mea calculated distribution matches exactly a theoretical 
distribution the chi-squared parameter is zero. 

Two methods can be used to calculate the chi-squared 
parameter for the goodness-of-fit test. In the first method, 
the €-axis of the PDF is divided into J intervals of equal 
size. The f,'s in this case represent the number of data 
points in each interval between divisions and f, varies 
across the distribution. In the second method, the PDF is 
g@ivided into J intervals of equal probability density, each 
containing approximately an equal number of data points in r.. 


Williams (1950) showed that to give the test sufficient 


mower the number of mimvervals, J, Shoulda De weiven by 
5 
peas 2(n-1) Cope 


10 





where n is the number of time samples. C, under Chie eau 
mapovhesis, is the cuvoifievelve for “the level of significance 


which as used in this analysis was taken as 0.05. 


yy 


VI. ANALYSIS OF DATA 


Hie analog data were digitized at 0.2 second intvervailc. 
The relationship between the upper limit of the bandwidth 


Genbe investigated, fa ewolely Ro lele” Felen io dialialcmmanlvaliviosmaes fe fei om she 
Af = ae. (Wo. ail,) 


Here r. is the Nyquist frequency which in this experiment 
was established at 2.5 Hz. This was sufficiently high to 
avoid aliasing or folding problems. The maximum lag time 
WeeecnOsen as consistent at 10 percent of the record engin. 
Maas results in twenty sabes of freedom for each. computed 
spectra. The 80 percent confidence limits for twenty 
degrees of freedom are between 0.71 and 1.60 of the measured 
spectral estimates according to the chi-squared distribution. 
Because the record lengths were different for each run, 
there was a different frequency resolution for each run. 

ive sirequency resolution for each. run is shown on the 


feepeclLive figures. 


A. VERTICAL VELOCITY MEASUREMENTS 

The predicted and measured vertical velocity energy- 
Gemcity spectra, and the phase angie y"ama Ccoherence=epectra 
of the wave height and vertical velocity were computed fOr 


runs two through six at the depths shown in Tabie [. 
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1. Predicted and Measured Vertical Velocity 
BMercy-Density Specura 


The vertical velocities were predicted using (4.4). 


The predicted vertical velocity-spectrum was determined from: 


jwce)|* = [oa Sink wiht)? yey]? (6.2) 
mae cerm in brackets is the square of the transfer functwon, 
(4.6), relating the vertical velocity to n. 

Measured vertical velocity energy-—density spectra, 
for Runs Two and Six, are shown in Figures (11) and (12). 

The principle energy band was from 0.04 to 0.37 Hz as deter- 
mined from the wave spectra. For Run Six, depth 5.70 meters, 
the theoretical values of the vertical velocity spectrum 
were within three to five percent of the measured vertical 
velocity spectrum values over a frequency range of 0.04 

to 0.34 Hz. The true error of the velocities is the square- 
root of the spectral differences and was approximately two 
percent. This is within the combined error of the wave and 
velocity meters; hence most of the measured energy-density 
within this bandwidth was wave induced with relatively 
little spectral energy-density being due to turbulence. 

On the other hand, for Run Number Two (the run of 
deepest depth at 14.62 meters) the absolute difference 
between the measured and theoretical spectra was small over 
the much narrower bandwidth of 0.05 to 0.19 Hz. The large 


absolute difference between the measured and theoretical 
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Figure ll. Measured and Theoretical Vertical Velocity 
Spectra for Run Number Two. Depth 14.62 Meters. 
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spectra at high frequencies for both runs indicates that 
beyond 0.19 Hz for Run Number Two and 0.34 Hz for Run 
Number Six, that most of the energy was due to turbulence. 

The decrease of wave induced motion with increased 
Gepth can also be seen by comparing the vertical velocity 
spectra for the two runs. For Run Number Six there was a 
peak in energy due to a three second wave component. This 
Peak was not present in the spectra for Run Number wo. 
That the energy present at this frequency and depth was 
Gue to turbulence can be seen by the large deviation 
between the measured and theoretical spectra. 

Like the spectral energy-density in the wave induced 
mere Of the spectrum, the turbulent enerey—-density lever 
was one order of magnitude less for Run Number Two compared 
Memrun Number Six. This indicates that turbulent motion 
@ecreased at the same rate as WavewanGduceda motion under the 
merLicular conditions of moderate swell, low winds, and 
intermediate water depth and that the intensity of the 
turbulence is associated with the rms value of the 


wave-induced motion. 

A measure of the total energy available in the 
measured spectra is given by the variance. The value of 
the variance, and therefore the total kinetic wave-induced 
energy present, decreased with increasing depth. There is 
almost an order of magnitude difference between the two 


Wellies Indiecavane G@arapid decrease inv enerey wath cep. 


as predicted by theory. 
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ce Coherence of the Wave Height and Measured 
Vertical Velocity 


ihe COonerence vepeerera. are also show? ine! eures 





(11) and (12). In general, the coherence was 0.8 or 
Beeaver throughout the range of signiiticant wave induced 
energy-density. This means that there is a hie wieceorreila— 
tion between the waves and wave-induced vertical motion and 
Meas relatively little of the energy in the Sierra wean 
emerey-density band was Que to turbulence, but dnostead. 

was wave-induced. 


3. Phase Angle Between Wave Height and the 
Vertical Velocity 


The phase angle theoretically leads wave hetehu by 
90 degrees. The two computed phase spectra in Figures (11) 
and (12) further reflect the accuracy of linear theory at 
both shallow and deep depths. 

Bordy (1972) encountered the greatest departure from 
linear wave theory in the computation of the phase angle. 
Much effort was expended during the course of this experi- 
Meme tO ensure that phase angle errors Gue to experimental 
technique were minimized. The Re BIOENG Ol tne weonseanec 
phase angle with depth rumen SUPPOreS Linear) ticory ele 
is felt that keeping the wave sensor directly over the flow- 
meter was the sole reason for the difference in the results 
obtained in this research as compared with those of Bordy. 


4, Comparison of Phase Angle and Coherence for 
Different Depths 


As can be seen from Figure (13), there is a definite 


relation between the phase angle for a given frequency and 
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Figure 13. Vertical Phase Spectra for Runs Two Through Six. 
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the depth of measurement. The 90 degree phase shift was 
observed to a period of 0.4 Hz at the shallowest depth, 
but only as far as 0.25 Hz for the deepest case. 

the upper frequency limit of high coherence also 
Meereases With ineréasing depth as illustrated in Figure 
(14). Both of these phenomena can be explained by linear 
mmcory. bamear theory. predicts decreasing particle motion 
Weel increasing depth due to @ passing wave and that bor a4 
Piven depth the lower frequency waves induce greater parti- 
ele motion. This is reflected in the coherence and phase 
munomws. The highest frequency @u whieh the Coherence 
approached one and at which the phase compared well with 
theory for all analyzed runs corresponded with the onset 


eresipnificant turbulent energy contribuyion. 


B. HORIZONTAL VELOCITY MEASUREMENTS 

mhe oredicted and measured WoC Dott Za tl aeavelocMl, 
energy-density spectra were computed for Run Number One. 
Hie ohase angle and coherence spectra were cCompuceq using 
the wave height and v-component of velocity. 

1. Predicted and Measured Horizontal Velocity Spectra 

The total theoretical horizontal velocity was 

computed from (4.3). The theoretical energy-density 
spectrum of the total wave-induced horizontal velocity was 


calculated from the expression 


aD? = Ue Ste) Ince? (6.3) 
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Figure 14. Vertical Coherence Spectra for Runs Two 
Oiveuch oda. 
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The term in brackets is (4.5), the transfer function 
melaving the total horizonval velocity to mn. Fhe total 
spectrum was found by calculating the individual u- anda 
V-component spectra and summing over the frequency band 
emeinverest. The spectra are iallustrated in semi-log formar 
tm Pigure (15). 

The difference hetween the spectra at the peak was 
found to be about 12 percent, giving an absolute value of 
meme difference of about 3.5 permeent in the Signiticant 


energy-density frequency band. 


a Coherence of Wave Hemeines and NMesaumed Toval 


Horazenvtal Vetoerr, 

The coherence (of n and the v-component of velocity) 
was calculated to be greater than 0.9, Figure (15), through- 
eue the significant energy frequency band of 0-03 to 020, 
me, this high coherence was due to tUhewiaet that the v— 
component of velocity was primarily wave-induced and the 
turoulent contribution to the spectrum was emall. 


3. Phase Angle Between Wave Height and the Total 
Horizontal Velocity 


The theoretical phase anple i ox Eye wave height 
and the v-component of velo 18 zero degrees. This was 
found to be true in this case as illustrated again in 
foure (15). 
©, PROBABILITY DENSITY FUNCTIONS AND VST Amis dGs- OF THE 

DATA SETS. 
A probability density function was computed for each 


Measured =dUantity Mer thetsi< runs.) 2nempom vs. Como tein 
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Figure 15. Measured and Theoretical Horizomeal Velocity 
Spectra for Run Number One. Depth 7.04 Meters. 
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a histogram of the data sets were plotted against both a 
Gaussian and a Crane cnenintes distribution. “Ihe eni-scquared 
goodness-of-fit test was computed for both theoretical 
@icoeribucvions. The statistics, variance, standard deviation, 
Seewness, and kurtosis were also computed for ali eighteen 
Meares celts. These resulesserecmmsvedsin table aie 
i. errobability Density sumer ons 

A typical PDF is @sheneman Pieure (to) meres c 
known that the ocean's surface has approximately a Gaussian 
eee rrabutIOn. However, Unem@mon-linearivies of (the suriace 
ierimoduce values of the Skewness anger kKUrvesie ther, deviave 
fee Gaussian and result inva Gram—-Charlier Gdistribuvion. 
toes 1S explained by the pnysiteal tact that waves tend vo 
werm short, sharp erests and long troughs. For the waves 
is Produces a positive skewness. indicavine more Values 
in the time series below mean water level than above. 

It follows that the waver indvced=partietve motaon 
should follow the same behavior. However, of the eighteen 
Gieaea Sets analyzed, three did not have positive skewness. 


These three were the velocity components u and w from Record. 
meve and w-Tfrom Record Tnrec. why  vntesisuso. despite tne 
mact, thal as3 the wave recordsr were Gpecir nem skeved, 

remains unexplained. A possible explanation could be that 


Bhis peculiarity is causea by internal waves which were 


prevalent during the experiment. 
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2. Goodness-of-Fit Test 

ive chi-squared goodness-of-fit test was used to 
Sompare the compuctec PDR to both the Gaussian and Gram— 
Charlier distributions for all eighteen data sets. 
Qualitatively, based on the smaller chi~squared values 
obtained from a comparison of the PDF with the Gram- 
Senor OVStri bution, 1G Can be svaved thavealimciehveen 
mecords more closely fit the Gram—Charliern qistripucvion 
than the Gaussian. The chi-squared parameter for all 


emeamveen Gata “sets is alsosiisetved in Table iy. 
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VIL. CONCLUSIONS 


The measurements in this research were made under condi- 
tions of moderate swell and low wind velocity. Wave height 
was measured with a penetrating wave staff and water 


pParcacle motion®measured Wath an eWeerromagnetie flowmeter. 


A. “@ENERGY—-DENSITY SPECTRA 

The wave-induced vertical velocity energy-density 
spectra were computed for six depths. The frequency of the 
eax Cnerpy-density remained, ceomeran’ with Gepcen. however: 
ticmovenitacant enercey=denscamearol! The wave—-imeuced mMovnon ar 
Nigh frequency decreased with increasing depth according to. 
Wace attenuation aspects of linear theory. © fhe magnavude of 
be enerey=-dcensity spectra also decreteecay vie 1ncreassne 
Sepun accordame to theory. It was moucdsuiay the lowes. 
Prequency at which turbulence comtrmeured Sienifteantly to 
wne spectra shifted to lower frequencies with increased | 
agepth. The intensity of the turbulent contribution decreased 
with depth and appears associated with the Vaetanee OF © ble 


wave-induced velocities. 


B. COHERENCE 

Coe of the wave height and wave-induced velocities 
was computed to be greater than 0.85 throughout the range 
Of significant enersy-densicy. Thewipper ereateney, Want or 
high comerence decreased with iAnereasiue depen. ine velve 


of the upper limit was established at the frequency where 


a 





turbulence became an important factor in the energy-density. 
The near unity values of coherence over the band of signifi- 
cant spectral-density imply a near linear process between 
Weve height and induced water particle motion. It is felt 
maae che small irregularities in coherence im the Significant 


emerey—density band were due to turbulence. 


C. PHASE ANGLE 

ine computed phase angle spectra of the wave heignt and 
induced water particle motions were consistent with the 
mmeory at all depths for both vertical and total horizontal 
velocity cases. This is the part of the analysis where 
Bordy (1972) encountered the greatest deviation from theory. 
It is imperative that the wave staff be directly over the 
mloymeter in order that there be no phase error. <A Similar 
decreasing upper frequency limit with increasing depth for 
consistent phase angle values was present as it was for both 


the coherence and the energy-density spectra. 


D. STATISTICAL INTERFERENCES 

The probability density functions of all measured 
quantities were qualitatively found to be better described 
py the Gram-Charlier distribution than by the Gaussian using 
the chi-squared goodness-of-fit test. The variance of the 
vertical velocities, a measure of the kinetic energy present, 
decreased with increasing depth which is consistent with 
linear theory. The skewness of the wave height frequency 


G2SturiburLron was found Go be positive in all cases tas <a 
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Consequence Of UnemiacteunareWwavesecena to Mave short. 
Sharp crests and elongated troughs. The values of skewness, 
With Tew exceptions, were found to be pesitiverior the waye— 
induced quantities as well. These exceptions might be 


explained by the presence of internél waves. 
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